Pulse Velocity in a Granular Chain 
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We discuss the applicability of two very different analytic approaches to the study of pulse prop- 
agation in a chain of particles interacting via a Hertz potential, namely, a continuum model and a 
binary collision approximation. While both methods capture some qualitative features equally well, 
the first is quantitatively good for softer potentials and the latter is better for harder potentials. 

PACS numbers: 45.70.-n,05.45.-a,45.05.+x 



The physics of a chain of particles interacting via 
a granular potential, i.e. a potential that is repul- 
sive under loading and zero otherwise, remains a chal- 
lenge despite a great deal of recent work on the sub- 

ject jaaiaaliaaaiiapEippiia. The- 

oretical studies of pulse dynamics in frictionless chains 
have relied primarily on numerical solution of the equa- 
tions of motion 0, |E H 13 13 13 E3 • Analytic work has 
relied on two rather different approximations with very 
little direct comparison between them. One approach 
is based on continuum approximations to the equations 
equations of motion followed by exact or ap pro ximate so- 
lutions of these approximate equations H, [Til [Jj| . This 
approach is expected to give useful results when the pulse 
is not too narrow, i.e., when the velocity distribution of 
the grains in the chain at any instant of time is rather 
smooth. The other approach is based on phenomenology 
about properties of pairwise (or at times three-body) col- 
lisions together with the assumption that pulses are suffi- 
ciently narrow to involve only two or three grains at any 
one time 0, 0, 0] . Among the interesting quantities 
one aims to calculate with these approaches is the pulse 
velocity. In turn, successful calculation of the pulse ve- 
locity requires a good understanding of the pulse width. 

The standard generic model potential between 
monodisperse elastic granules that re pel upo n overlap ac- 
cording to the Hertz law is given by [2(J, |21| 

V(8 k>k+1 ) =mt k+ v S<0, 

(1) 

V{5 k ,k+i) =0, 6>0. 

Here 

£fe,fe+l = Vk - Vk+l, (2) 

a is a constant that depends on the Young's modulus and 
Poisson's ratio, and yk is the displacement of granule k 
from its equilibrium position. The exponent n is 5/2 for 
spheres, it is 2 for cylinders, and in general depends on 
geometry. 

The displacement of the fc-th granule (fc = 1, 2, . . . , L) 
in the chain from its equilibrium position in a frictionless 



medium is thus governed by the equation of motion 

+a(y fe _! - ykT^eivk-i - Vk). (3) 

Here 9(y) is the Heavyside function, 9(y) = 1 for y > 0, 
9{y) = for y < 0, and 9(0) = 1/2. It ensures that the 
particles interact only when in contact. Note that for a 
finite open chain the first term on the right hand side of 
this equation is absent for the last granule and the second 
term is absent for the first. 

Initially the granules are placed along a line so that 
they just touch their neighbors in their equilibrium posi- 
tions (no precompression), and all but the leftmost parti- 
cle are at rest. The initial velocity of the leftmost particle 
is vq (the impulse). In terms of the rescaled variables 

(2 \ 1/™ i / 2 \ i/ n 

— x k , T=— — M t, (4) 
a J v \ a J 

Eq. can be rewritten as 

x k = -{x k - x k +i) n ~ 1 9(x k - x k +i) 

+(£fc-i - x k ) n ~ 1 9(x k ~i - x k ), (5) 

where a dot denotes a derivative with respect to t. In the 
rescaled variables the initial conditions become x k (0) = 
Xk (0) = 0, Vfc ^ 1, xi(0) = 0, and ±i(0) = 1. 

When n > 2 an initial impulse settles into a pulse that 
is increasingly narrow with increasing n, and propagates 
at a velocity that is essentially constant and determined 
by n and by the amplitude of the pulse. For n = 2 the 
pulse spreads in time and travels at a constant velocity 
independent of pulse amplitude. In the latter case there is 
considerable backscattering that leads to backward mo- 
tion of all the granules behind the pulse, whereas the 
backscattering is minimal for n > 2 |, [z| • The pulse 
is a completely conservative solitary wave in the limit 

71 — » OO. 

Three features determine pulse dynamics in these 
chains: 

1. The power n in the potential; 

2. The absence of a restoring force; and 
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FIG. 1: Pulse velocity as a function of the power of the po- 
tential. The stars represents the numerical simulation results, 
the continuous line the binary collision approximation, and 
the broken line the continuum approximation. 



FIG. 2: Relative error (CV — Cs)/Ct in the pulse velocity 
as a function of the power of the potential. Here, the the- 
oretical pulse velocity Ct, is either the one obtained from 
the continuum (circles) or from the binary collision (squares) 
approximation. Cs is the pulse velocity obtained from the 
numerical simulations of the chain. 



3. The discreteness of the system. 

Recently, we discussed the role each of these features in 
the continuum approximation and extended previous re- 
sults to include viscosity [22] ■ Not only does this approx- 
imation work very well for the n = 2 case H, l2^ j. but, 
corroborating Nesterenko's theory 0, we found that the 
continuum approximation works surprisingly well for the 
prototypical grain, namely spheres. However, a discus- 
sion of the reasons for this agreement seems to be lack- 
ing. On the other hand, models based on binary interac- 
tions have also been proposed in order to study a chain 
of spherical and other grains. Wu's independent-collision 
model [LH focuses on a chain of tapered grains. From en- 
ergy and momentum conservation considerations, work- 
ing in the n — > 00 limit, he shows that his simple model 
captures the qualitative behavior observed in simulations 
for spheres. Subsequently this model was phcnomeno- 
logically extended and compared with experimental re- 
sults [l|. 

Herein we address the question of the applicability of 
both the continuum approximation and a binary inter- 
action model through the analysis of the velocity of the 
signal propagation as a function of the power of the po- 
tential. For the former case, the pulse velocity, C c (n), 
can be written as [22j 



C c (n) 



n 2 (n - 2) 



2(n + 2) v /^/(^ 



(6) 



where 



1(1) 



2 *r 2 (¥) 



r(z + i) 



(J) 



On the other hand, for the binary collision approxi- 
mation the set of equations (J5J reduces to two coupled 
equations which may be decoupled by defining the nor- 
mal mode variables z± = x\ ± X2- In particular, we have 



= -2z n ~ 



(8) 



This is precisely the equation of motion for one parti- 
cle subjected to a potential V(z) = 2z n /n. Furthermore, 
the initial conditions for the original variables lead to 
z_(0) = and i-(0) = 1. Hence, from energy conserva- 
tion we have: 



1 1 9 

-z 2 (t)-- = --z n (t). 



(9) 



Consequently, when the two particles have the same ve- 
locity (i_ = 0), their compression will be a maximum 

and given by z m — (j) 1 ^™. Therefore, the time neces- 
sary for the "pulse" to travel from the first to the second 
particle is 



T h (n) 



yj\ - \z n _jn 



zdz_ 



(10) 



Explicitly integrating Eq. (|10|) . we can write the pulse 
velocity C^in) — 1/T b (n), as 



C b (n) 



,1/2 



l/n 



r(i + i) 



(ii) 



Our comparison is thus between Eqs. © and 

Qualitatively both approximations give the same re- 
sult: for n > 2 the pulse velocity decreases with n, at- 
taining its minimum value for n~5, and then increasing 
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and saturating for large n (see Fig. ^| . Quantitatively, 
however, they differ appreciably. For instance, for large 
n, Cb — > 1 while C c — > 0.883 .... In Fig. [3 we present the 
relative error of our numerical simulation of the chain 
as compared with each theory. From this figure, it is 
easy to see that while the continuum approximation gives 
very good results for small n, its predictions are poor for 
n > 3.5. On the other hand, the binary collision approxi- 
mation is extremely accurate for n > 3.5 but not accurate 
for small n. The quantitative agreement of each of the 
two analytic results at their respective n extremes with 
the numerical simulations is seen to be excellent. Again, 
we stress that these results are a reflection of the behav- 
ior of the pulse width. At smaller n the pulse is relati vely 
broad, the velocity pulse covers a number of grains [22|, 
and a continuum approximation captures the pulse con- 
figuration and speed very well. At larger n the pulse 
becomes very narrow, discreteness effects dominate the 



behavior, and an approximation that assumes that two- 
grain collisions dominate the pulse behavior reproduces 
the pulse velocity extremely accurately. 

In summary, Nesterenko's continuum approximation 
gives quantitatively accurate results for the pulse velocity 
for relatively soft potentials, n < 3.5 (which includes the 
generic cases of cylindrical and spherical grains), while 
the binary collision model is quantitatively correct for 
relatively hard potentials, n > 3.5. 
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